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Abstract. This paper is the first in a series devoted to the hard X-ray whole sky survey performed by the 
INTEGRAL observatory over seven years. Here we present an improved method for image reconstruction with 
the IBIS coded mask telescope. The main improvements are related to the suppression of systematic effects 
which strongly limit sensitivity in the region of the Galactic Plane (CP), especially in the crowded field of the 
Galactic Center (GC). We extended the IBIS/ISGRI background model to take into account the Galactic Ridge 
X-ray Emission (GRXE). To suppress residual systematic artifacts on a reconstructed sky image we applied 
nonparametric sky image filtering based on wavelet decomposition. The implemented modifications of the sky 
reconstruction method decrease the systematic noise in the ^ 20 Ms deep field of GC by ^ 44%, and practically 
remove it from the high-latitude sky images. New observational data sets, along with an improved reconstruction 
algorithm, allow us to conduct the hard X-ray survey with the best currently available minimal sensitivity 3.7 x 
10~^^ erg s~^ cm~^ ^0.26 mCrab in the 17-60 keV band at a 5cr detection level. The survey covers 90% of the 
sky down to the flux limit of 6.2 x 10~^^ erg s~^ cm~^ (~ 4.32 mCrab) and 10% of the sky area down to the flux 
limit of 8.6 x 10~^^ erg s"^ cm"^ (- 0.60 mCrab). 

Key words. Surveys - X-rays: general - Galaxy: general - Methods: data analysis - Methods: observational - 
Techniques: image processing 



1. Introduction 

Since it s launch in October 2002, the INTEGRAL obser- 
vatory ( Winkler et al. |2QQ3 ) has gathered a huge observa- 
tional data set allowing us to perform the most sensitive 
hard X-ray survey to date. The main scientific results and 
source catalogues have been reported in many relevant pa- 



pers concerning partial sky coverage (e.g. Revnivtsev et al. 



2003; Molkov et al. 2004; Krivonos et al. 2005; Revnivtsev 
et al. 2006a, Bird et al., 2004 



2006 



et al. 



,2006, 2007, Bassani et al. , 
Bazzano et al.||20Q6 ) and full sky studies (Krivonos 



..2007b| ISazonov eTari[2QQ7| [Beckmann et al.||2009[ 
Bird et al.| |2010l). 

Recently, great progress in surveying the hard X-ray 
sky was achieved with the Burst Alert Telescope (BAT; 



* Based on observations with INTEGRAL, an ESA project 
with the instruments and science data centre funded by 
ESA member states (especially the PI countries: Denmark, 
France, Germany, Italy, Switzerland, Spain), Czech Republic 
and Poland, and with the participation of Russia and the USA. 



Barthelmy et al.|[2005| ) at the Swift observatory (Gehrels 
et al.||2004 |).The Swift /BAT survey provides very homo- 



geneous sky coverage in the 15-195 keV energy band with 
a current maximum sensitivity of 2.2 x 10~^^erg s~^ cm~^. 
The distribution of survey sensitivity peaks in the ex- 
tragalactic sky. The survey results and source catalogues 
have been reported in papers bylTueller et al. (2009) and 



Cusumano et al. (2010). As seen from the large sample of 
detected Active Galactic Nuclei (AGNs), the results of the 
Swift /BAT survey are very valuable for extragalactic stud- 
ies. However, due to the relatively poor angular resolution 
of the instrument, its capabilities in the Galactic plane 
and especially in the Galactic Center regions are limited. 
On the other hand, the sky coverage by the Swift /BAT 
survey is nearly uniform, therefore only a small fraction 
of its total operational time was devoted to observations 
of the Galaxy. 

Contrary to Swift, the INTEGRAL observatory pro- 
vides an all-sky survey with exposure more concentrated 
in the Galactic Plane, having a typical limiting flux of less 
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than 1.43 x 10~^^erg cm~^ (1 mCrab) in the working 
energy range 17-60 keV. With an angular resolution al- 
most twice as good as Swift /BAT, one can effectively dis- 
entangle sources in such crowded regions as the Galactic 
Center. This makes the Swift /BAT and INTEGRAL sur- 
veys complementary to each other. 

INTEGRAL has already accumulated a lot of expo- 
sure time in the direction of the Galactic Plane with a 
maximum of ~20 Ms of nominal time in the direction of 
the Galactic Center. However, the growing exposure time 
devoted to the Plane and the Center of the Galaxy is not 
reflected by a corresponding increase in survey sensitivity. 
Observations in these regions are strongly affected by the 
systematics related to the crowded field of the GC and 
strong Galactic X-ray background radiation. 

In this work we address the question of improving 
the sensitivity of the ongoing INTEGRAL hard X-ray 
survey. In section Q we discuss several aspects of the 
sky reconstruction method of the IBIS coded-mask tele- 
scope ( Ubertini et al. [2003 ). In section ([3| we implement 
Galactic background corrections to the sky reconstruction 
method. Section Q introduces a modified sky reconstruc- 
tion method with sky image filtering procedure based on 
d trous wavelet decomposition. The properties of the re- 
sulting all-sky survey are presented in section ([5|. 

Throughout the article the exposure will be expressed 
taking into account instrumental dead time, i.e. showing 
the effective exposure time, rather than the total exposure. 

2. General sky reconstruction method 

IBIS is a coded aperture imaging telescope. The sky is 
projected on to the detector plane through the transpar- 
ent and opaque elements of the mask mounted above the 
detector plane. Generally, the sky reconstruction is based 
on the deconvolution of the detector image with a known 
mask pattern. We (EC) developed IBIS/ISGRI sky recon- 
struction method and partially described it in our previous 
pubhcations ( Revnivtsev et al.|2QQ4| Krivonos et al.|2QQ5 



2QQ7a^b}. T he basic idea we used is presented in |Fenimore 
Cannon | ( |1981 ) and Skinner et al. (1987). For the stan- 
dard IBIS/ISGRI analysis we refer reader to the paper by 
Goldwurm et al.| ( 2QQ3J . Here we outline only those steps 
that are essential for the present study. 

2.1. IBIS telescope coding aperture 

The quality of the reconstructed sky image directly de- 
pends on our understanding of the coding procedure. For 
example, the mask supporting structure can significantly 
reduce the number of low energy photons passed though 
the open mask elements. Other assembly elements like 
screws, plates and glue strips attaching the IBIS mask 
to the supporting structure also block the incoming pho- 
tons modifying the shadow cast by a point source on the 
detector. 

In our package we have implemented the best known 
configuration of all the known elements of the telescope. 



Fig. 1. Sky image of a 2.7 ks observation of the bright 
source Crab Nebula. The source is in the center of the 
field of view. The periodic mask pattern leads to many 
false sources on the reconstructed sky. The most signifi- 
cant false peaks are located in vertexes of a ~ 10.5° square 
around the source, oriented in detector coordinates as il- 
lustrated by the yellow square. The image is shown in 
significance with a squared root color map ranging from 
to 25. Such color schemes are used for all the sky im- 
ages in this paper in order to emphasize sky background 
variations. The black and blue colors correspond to pixel 
values from to 2. The red pixels have values of around 
5. The yellow to white color transition corresponds to 15 
and more. Fig. [2] demonstrates an image profile extracted 
from the 1° x 29° green region. 




FOV offset, degrees 

Fig. 2. Image profile extracted from the green region 
in Fig. [l] The profile demonstrates the response of the 
IBIS/ISGRI imaging system to the point source. 
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which make contributions to the shadowgram. We could 
not find the exact size of the assembly screws which at- 
tach the mask, plates and glue strips. The approximate 
parameters of these elements were found by comparing 
real and model detector shadowgrams illuminated by a 
strong source at different azimuthal angles. 

2.2. Detector exposure 

In general the detector image, produced when observing 
M point sources with an underlying fiat sky background, 
can be represented by the superposition of shadow pat- 
terns of sources ("pixel illumination fraction", PIF) and 
the detector background map B: 



M 



D = ^f,PIFi^kBB, 



(1) 



i=0 



where fi is the source fiux, B = Bex b + ^det represents a 
detector background map, containing photon counts from 
Cosmic X-ray Background (CXB), and detector instru- 
mental noise. We assume, that the detector illumination 
by CXB and intrinsic ISGRI background have a similar 
pattern, and can be merged into the background map B. 
The last is estimated by an accumulation detector im- 
age over a large number of observations without strong 
sources in the field of view. Obviously, background map 
B should contain the current detector background pattern 
due to the long-term variation of background environment 
related to the Sun and Cosmic rays (Lebrun 2005). For 



a given observation, we use a background map accumu- 
lated during the nearest set of extragalactic observations. 
Typically, we construct a new background map for every 
- 50 - 70 spacecraft orbits (150-200 days). 

In the case when two sources located close to each 
other (at separation comparable with IBIS/ISGRI angular 
resolution) , the direct solution of Eq. [l] for fiux fi can be 
unstable giving results with infinite errors. In other words, 
the detector count rate in a given sky direction can be 
explained by two (or more) sources having any absolute 
fiux. 



2.3. Replicated mask pattern 



IBIS mask has replicated patterns (see Reglero et al.|2001 



Ubertini et arf2003, .Goldwurm et al.|2003D . This pattern 



of the mask has an advantage because, ideally, it has a 
much narrower point spread function and fiat side lobes 
in the central part of the reconstructed image (Fenimore^ 
fc Canno'n~||1978 ). But at the same time it causes serious 
complications due to the presence of the very significant 
side peaks of the point spread function. This means that 
a simple deconvolution algorithm sees "ghost" sources at 
certain sky positions (see Fig. [l] and [2]), related to the 
position of the real source and the size of the replicated 
mask pattern. 

The sources in a variety of sky positions within a field 
of view create shadows with similar patterns, which causes 



uncertainty of source fiux determination (the direct solu- 
tion of Eq. [l]is impossible). Unfortunately, this situation 
is not rare in the crowded field of the Galactic Center, as 
shown in Fig. [3j 

2.4. Iterative removal of cataloged sources 

In order to obtain a good quality all- sky map, suitable to 
searching for new weak sources, the ghosts of known bright 
sources have to be removed. This is done during the re- 
construction of images of individual observations {ScWs) 
with an effective exposure of 1 — 3 ks. Instead of a blind 
search for bright sources in each individual observation we 
use a catalogue of known sources to control the removal 
of ghosts. Indeed, for many regions of sky the final map 
is the result of stacking hundreds and thousands of indi- 
vidual observations. Therefore, a relatively weak object, 
far too faint to be detected in an individual observation, 
may appear as a very significant source in the final map. 
Since the amplitude of the ghosts scales with the intensity 
of the true source, it is clear that ghost removal should be 
applied even to objects which are too faint to be detected 
in individual science windows. For crowded fields (like the 
Galactic Center region) this implies that ghosts of some 
100 sources should be removed in individual observations. 

The whole procedure requires not only the list of 
sources to be removed from the detector image, but also 
a sequence of Iterative Removal Of Sources (IROS, see 
Goldwurm et al.| ( |2003D ; [Krivonos et^ ( |2005D ). It is ex- 
pected (and confirmed by direct tests) that the brightest 
(most significant) objects have to be removed first, since 
the source fiux / is evaluated assuming that there is only 
one source in the field of view. The significance of the 
source detection is evaluated by reconstructing an image 
prior to iterative source removal and checking the fiuxes 
at the positions of cataloged sources. The list of objects 
ranked according to their significance is then used as in- 
put for iterative source removal procedure. This poses the 
problem of ranking weak sources, since their fiux (and 
ranking) is determined with a large uncertainty in the in- 
dividual science window (e.g. the fiux from the source can 
be negative). We made several tests with various ranking 
schemes for weak (less than 3a detection) sources, check- 
ing the RMS of the final maps and fiuxes of cataloged 
sources. The final scheme implemented in our analysis uses 
the absolute value of the source detection significance to 
rank the order of source removal. 



2.5. ISGRI pixel filtering 



As described in Krivonos et al. (2007b) the hot and dead 



ISGRI detector pixels were screened from the analysis. 
This was done using quite crude filtering criteria and some 
noisy pixels may still be present on the detector shadow- 
gram. They may not be visible on an individual detector 
image, but can be revealed by those accumulated over sev- 
eral observations. 
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Fig. 3. INTEGRAL/IBIS hard X-ray (17 — 60 keV) map of the sky region around the Galactic Center. The green 
squares demonstrate the relative sky positions of false side peaks ("ghosts", see Figjl]) of selected bright sources. Due 
to observational constrains IBIS FOV was mainly co-aligned with axes of the equatorial coordinate system ("N-E" 
notation). As a result, there are a number of bright sources which permanently appeared in mutual "ghost" positions. 
This fact leads to an additional uncertainty of the source flux determination. The following source pairs can be affected 
by mutual flux interplay: GX 340+0 and XTE J1743-363, XTE J1701-462 and GX 349+2, GX 9+9 and GX 13+1, 
and the triplet GX 1+4, XTE J1818-245 and GX 17+2. 



To estimate the effect induced by "noisy" pixels, we 
simulated a detector image for a typical ScW exposure of 
2 ks, and inserted one pixel in an arbitrary position ex- 
ceeding the mean detector count rate (40 cnts/pix/Scl^) 
by a factor of ~ 2.5 (^ 9.5 standard deviations). For a typ- 
ical ScW such a weak pecular pixel introduced negligible 
systematic noise and the reconstructed sky was dominated 
by Poisson statistics. However, with increasing exposure, 
the effect became more significant. We accumulated the 
mosaic image of a 280 ks staring observation of NGC 4151 
with a simulated detector and one hot pixel. When the 
position of a noisy pixel was randomly distributed on a 
detector in every staring observation, the total mosaic did 
not contain any significant systematic residuals. When a 
noisy pixel was fixed on the detector, the standard devia- 
tion of reconstructed sky was ~ 20% higher, than without 
a hot pixel on the detector. 

To perform additional ISGRI pixel cleaning, we fol- 
lowed the general approach also employed by Eckert et, 



al 



(2008). We stacked detector images obtained dur- 



pected that 99.7% of the pixels have a value in the range 
[-3cr,+3cr]. The ISGRI detector contains 128 x 128 pix- 
els which gives us ~ 50 of them with an expected value 
greater than 3a. However, we typically detected 200 — 300 
deviations from zero to larger than 3a. We removed these 
pixels from further analysis. Filtering done in this way re- 
duces the detector area by ^ 2%, and has a minor effect 
on fiux from serendipitous faint sources. 



3. Galactic background 

During observations of the Galactic Center region the IBIS 
field of view contains many discrete sources (Fig.[4|.But m 
contrast to high galactic latitude observations, the under- 
lying sky background is not flat. From early X-ray obser- 
vations we know, that the Galaxy reveals itself as a strong 



diffuse emitter (e.g. Worrall et al. 1982). The morphol- 



ogy of the Galactic X-ray background at energies above 
20 keV is now relatively well known. As shown in recent 
RXTE and INTEGRAL investigations (iRevnivtsev et al. 



ing the spacecraft orbit after removing flux from known 2006b Krivonos et al. 2007a), the X-ray background is 



X-ray sources (Sect. 2.4), detector background map B 
(Sect. [2^2] ), and Galactic X-ray background (Sectj3|. The 
distribution of pixels on a stacked detector image was 
well described by Gaussian with zero mean. Thus we ex- 



traced by the near infrared brightness of the Galaxy (blue 
contours in Fig. |4|. We will refer to the Galactic X-Ray 
Background later as "Ridge" emission or GRXE (Galactic 
Ridge X-ray Emission). 
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Fig. 4. Two relative alignments of GP in the IBIS FOV, 
left centered and the right shifted along the galactic lat- 
itude. The blue contours are isophotes of the 4.9/im sur- 
face brightness of the Galaxy (COBE/DIRBE) revealing 
the bulge/disk structure of the inner Galaxy. The NIR 
brightness of the galaxy traces the hard X-ray Ridge emis- 
sion. The small and large rounded squares on each plot 
demonstrate the full and partial coded areas respectively. 
The points show sky positions of the hard X-ray sources 
detected on the 20 Ms time-averaged map (Fig. 11|3 ). 
A signal-to-noise ratio of the sources is shown by black 
(5 - 10), blue (10 - 30), and red (> 30). The difference of 
detector illumination by the Ridge is shown in Fig. [5] 
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Fig. 5. The contribution of Galactic X-ray Background to 
the ISGRI detector count rate for different alignments of 
GP in the IBIS FOV (see also Fig.[4|. The detector count 
rate normalization is estimated from actual observation of 
GC in the energy range 17-60 keV. The corresponding im- 
ages m Fig. [4] and |5] are co- aligned in detector coordinates. 
The Ridge, appeared in the given FOV corner, illuminates 
the opposite corner in the detector image. 



The measured 17-60 keV GRXE intensity per IBIS 
FOV reaches 200 mCrab in the region of the Galactic 
bulge. Such strong emission will not appear on a decon- 
volved IBIS/ISGRI sky, because the coded mask tech- 
nique is not suitable for building an image of sources 
with angular sizes larger than the telescope's angular res- 
olution (12^). However, the presence of Galactic X-ray 
Background in the IBIS FOV strongly affects the over- 
all shape of the ISGRI detector shadowgram, which leads 
to appearence of systematic noise on the reconstructed sky 
images. 

The detector exposure by the point source and GRXE 
is different. The point source forms a shadow on the de- 
tector. Each exposed pixel contains an approximately con- 



stant source count rate, therefore the source shadowgram 
is flat (in illuminated pixels) . This is explicitly assumed in 
Eq.[l} When IBIS FOV contains the Galactic Ridge emis- 
sion, one detector pixel "sees" different parts of the Ridge 
through the mask's open elements. The resulting detector 
image is convolution of the Ridge sky brightness distribu- 
tion with the IBIS mask and the collimator response. 

In order to model the contribution of the Ridge com- 
ponent to the detector we have to assume some predefined 
sky surface brightness distribution. The Ridge emission at 
energies 17 — 60 keV is well traced by the near infrared 
brightness of the Galaxy. The volume density distribution 
has been intensively investigated by COBE/DIRBE (see 
e.g. pwek et al.|1995D. We use the Galactic disk and bulge 
model from |Revnivtsev et al. (2006b), which describes the 
COBE/DIRBE data in the simplest way. The model was 
taken from the near-infrared data and renormalized to fit 
X-ray intensity measured by INTEGRAL. We produced 
sky brightness maps by integrating the model volume den- 
sity through the line of sight. 

The contribution of the Ridge emission to the ISGRI 
detector can be implemented with the so called "gray 
mask" approach: the given detector pixel sees the Ridge 
emission through the open mask elements in a solid angle 
limited by the telescope collimator. 

To demonstrate Ridge detector exposure, we projected 
GRXE sky brightness onto the detector using two posi- 
tions relative to the IBIS FOV. In position "A" the cen- 
tral part of the Ridge emission is placed in FOV on-axis 
(spacecraft pointing / = 0°,6 = 0, Fig. [4j left), and posi- 
tion "B" when the Ridge is 7° away from the FOV center 
(/ = 0°,6 = 7°, Fig. right). The Ridge detector illu- 
mination implemented with the "gray mask" is shown in 

Figig 

As clearly seen from the modeled detector images, the 
Ridge emission introduces large scale variations on the 
ISGRI detector. When IBIS is centered on the Galactic 
plane (position "A"), the Ridge contribution to the de- 
tector image has a relatively flat shape with some cur- 
vature. In contrast detector variations have a significant 
gradient when the GP intersects FOV several degrees away 
from its center (position "B"). That configuration strongly 
warps the detector image, which obviously will affect the 
source flux estimates and, consequently, the accuracy of 
the source shadowgrams removal. 

The impact of the presence of the Ridge emission in the 
IBIS FOV to the reconstructed sky image is significant. 
The comprehensive demonstration of this effect directly 
on the GC data is not possible due to more serious effects 
related to the very complicated detector exposure to many 
bright sources. That is why, in order to show the Ridge 
contribution in the reconstructed sky, we used the rela- 
tively simple and clean 280 ks observation of NGC 4151 
for reference (Fig. [sj left). The artificial detector Ridge 
component with actual normalization measured in the GC 
(position "B") was added to the detector image of every 
spacecraft observation. Since we used staring observations, 
this operation was equivalent to placing Ridge on the sky 
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Fig. 6. Left: The sky mosaic of 280 ks staring observation 
of NGC 4151. The Ridge contribution (Fig. [sj right) was 
artificially added to the detector plane for each individual 
observation. For reference see the left image in Fig. [Sj 
Right: Signal-to-noise ratio distribution of a number of 
pixels in a reference and current sky images shown by red 
and blue histograms, respectively. The long dashed line 
represents the normal distribution with unit variance and 
zero mean. 



7° away from the NGC 4151. The quality of the final mo- 
saic (Fig. [6j left) was very poor and the signal-to-noise 
distribution of pixels had strong non-Gaussian wings (blue 
histogram in Fig.[6|. The standard deviation of image was 
1.7, in contrast to 1.3 of the referenced sky. Thus, we can 
conclude that the Ridge emission can introduce significant 
systematic noise to the reconstructed sky. However degree 
of image worsening depends on exposure and pattern of 
the observations. 

Let us consider the typical IBIS/ISGRI observation 
of the Galactic Center region (position "A"), as the most 
representative observation containing many bright sources 
and strong Ridge emission in the FOV. The detector image 
is shown in the left frame of Fig. [7[ 

Firstly, the catalogue of predefined source positions 
provides about a hundred objects in the FOV. Among 
them, 6 — 8 bright sources are usually detected with 
S/N > 5 in the individual observation of ~2.5 ks expo- 
sure time. By summing up all PIFs of these sources, one 
can show that practically all detector pixels are illumi- 
nated by at least one source (see Fig.[7) left). This means, 
that any measurements of source fiuxes in the GC will be 
affected by the correlations of source shadowgrams, espe- 
cially for those having replicated patterns due to the peri- 
odic mask elements. Furthemore, there are no source-free 
pixels to estimate the background count rate. Flux mea- 



surement, based on the balance matrix (Krivonos et al. 
2005 Goldwurm et al.||2003 ) is not accurate, because the 
background map for a given source is strongly affected by 
other sources. The general reconstruction algorithm works 
"as is" , trying to estimate the fiuxes of the brightest source 
and leaving many systematic artifacts on the detector, and 
consequently, on the resulting sky image. 

The situation is further complicated by the fact, that 
practically all ISGRI pixels are illuminated at the same 
time by the Ridge emission (see Fig. [5|. This leads to a 




Fig. 7. Left: The ISGRI detector image of the Galactic 
Center (/ = 3.9°, 6 = 2.1°). Right: The shadowgram of 
a sky region containing only one bright source Sco X-1 
(/ = 4.5°, 6 = 30.0°). The observations were separated 
by a 10 hour time interval when the background count 
rate had not significantly changed, i.e. the background 
map for both images is the same. In the right image the 
background and source counts are clearly distinguishable. 
This leads to the straightforward application of the stan- 
dard coded mask reconstruction algorithm (see Sect. [2|. 
The left image contains many overlaping shadowgrams of 
bright sources. The assumed background model is not ac- 
curate and the general method is confused, (see text) 



high correlation between shadowgrams of a manifold of 
sources and the Ridge. 

Summing up all mentioned effects, we conclude, that 
Eq. [l] based on simple background map B is not suitable 
for Galactic observations, and the iterative source removal 
is not valid when the detector map is warped by the Ridge. 
Generally, sky reconstruction is highly complicated in the 
case of Galactic Center observations. To a large extent, 
the INTEGRAL/IBIS/ISGRI survey sensitivity is limited 
by systematic uncertainties, and new exposures in the GC 
provide only minor effects on the total sensitivity. 

4. Modified sky reconstruction method 

The sky reconstruction method for IBIS/ISGRI should be 
able to split a detector image into the following layers: 

— source PIFs. Characterized by a fiat count rate and 
the pattern of the projected mask. 

— illumination by the Ridge. This layer has certain low- 
(spatial)-frequency variations. The characteristic vari- 
ations depend on the telescope orientation relative to 
the Galactic Plane and Center. 

— ISGRI background layer. This contains the CXB fiux 
and instrumental background, and has a smooth count 
rate. 

all these layers correlate in different degrees with each 
other due to the shared detector area. We modify Eq. [l] 
to accommodate the Ridge component R: 



M 



D = Y.fiP^^^^^BB^kRR 



(2) 



z=0 
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The direct simultaneous solution of Eq.|2]is practically 
impossible before the components are orthogonalized. At 
first, we tried to reduce the mutual correlation between 
layers using different spatial features of the layers and the 
short list of sources (to do this, we selected M5 bright 
sources with S/N > 7 which appeared in the FOV). We 
fit the relation ([2| with the modelled Ridge contribution 

source PIFs, and background map B. The fitting pro- 
cedure was unstable, producing inadequate estimates of 
the layer normalizations. We noticed high uncertainty in 
determination of Ridge and Background map components, 
respectively kR and /c^, due to its high mutual correlation. 

One possible way to further disentangle the Ridge layer 
from others is to use different spacecraft orientations rela- 
tive to GP. To do this, we attempted to simultaneously fit 
the nearest Q observations (ScWs) before and after the 
current one: 



Mb 



(3) 



i=0 



where index j is running in the range [— Q, +Q] relative 
to the current observation j = 0. M5 denotes the num- 
ber of bright sources which appeared in the FOV of every 
j G [— +Q] observation with S/N > 7. The on- axis di- 
rection of all selected observations is requested to be inside 
15° radius around the position of current observation. 

We fitted large data sets using Eq.[3]for different Q > 
and found the Ridge component kn much better con- 
strained with respect to the other layers. Guided by the 
stability of the fitting procedure, we chose Q = 6 as an rea- 
sonable selection criteria. Generally, the procedure is more 
stable for more spatially scattered observations in the di- 
rection perpendicular to the GP. Among the available 
INTEGRAL observational patterns we have found that 
the best pattern for this approach is the Galactic Center 



Deep Exposure (GCDE, Core Program, see e.g. Winkler et 
al. II999I) and the Galactic Latitude Scans (PI Sunyaev). 
The fitting procedure running on large data sets of the 
usual 5x5 pattern (see ISOC Newsletter #12, September 
2004) and the Galactic Plane Scans (GPS, Core Program, 
see e.g. Del Santo et al.|2003 Rodriguez et al.|2QQ3J is less 
stable, but still provides valuable results. 

We should note, that the contribution of the ith source 
{fiPIFij) in Eq. |3]is estimated under the assumption of 
the constant flux fi during the considered time interval. 
In a similar way, the detector background count rate ks 
is also considered constant. For the chosen value of Q = 
6, the maximum number of observations is 13, which is 
in total ^ 25 ks. Generally, for such a time interval, the 
majority of galactic X-ray sources do not vary by a factor 
of more than ~2, except during the outburst activity. 

Actually, the modified sky reconstruction algorithm 
described here is used only for constraining the Ridge com- 
ponent and its following subtraction. After this step, the 



IROS procedure (Sect. 2.4) is applied in the usual manner 
on the detector plane of the given observation (ScW). 



Note, that the employed procedure still doesn't allow 
us to completely resolve the problem of highly correlated 
shadowgrams on the detector in the case of the Galactic 
Center observations, but it at least reduces the systematic 
residuals on mosaic images introduced by Ridge emission. 



4.1. Removing systematic residuals from sky images 

After removing the source shadowgrams and the back- 
ground from the detector we still see the systematic ef- 
fects on sky images. These residuals can be clearly seen on 
deep extragalactic observations in so called staring mode. 
In these observations orientation and the roll angle of the 
telescope is fixed. The sky images are stacked (pixel-to- 
pixel) in detector coordinates. All systematic residuals, 
not visible on images of individual observations (with an 
exposure of ~2 ks) are amplified on the stacked image. The 
sky region around NGC 4151, was reconstructed with the 
help of the general deconvolution algorithm from the data 
collected in the staring mode, as shown in the left image 
of Fig. |8] The characteristic chessboard-like squares and 
ripples are clearly seen in the fully coded (central 10° x 10° 
square) and partically coded (outer parts) field of view. 

Due to the absence of bright sources in the field of view 
during these staring observations, it is clear that the pres- 
ence of the patterns seen in Fig. [8] does not depend on the 
accuracy of our model of source shadowgram. We cannot 
also attribute visible systematic artifacts to the detector 
noisy pixels, because they were filtered out (Sect. [2^ . 

The major origin of these patterns on the sky is the 
limited knowledge of the background pattern on the de- 
tector. It may consist of several unaccounted for parts, like 
the unexpected variations of pixel gains, the effective life- 
times, and lifetimes of the detector modules (the ISGRI 
detector has 8 modules of pixels, which in many cases 
change, see e.g. Fig. [5|. Judging from a particular pattern 
of residuals on the sky, the major effect is due to the inac- 
curate estimate of efficiency and the energy band-passes 
of the detector pixels. 

In the context of the general sky reconstruction 
method, it becomes extremely difficult to correct this 
problem due to the continuous degradation of detector 
pixels and the variation of the background environment. 
Therefore, while understanding the origin of the residual 
structures on the sky we decided to implement an alter- 
nate solution to this problem. 

As the spatial scale of systematic artifacts on the sky 
is significantly different from that of the point sources, 
we have implemented the wavelet-based image filtering 
procedure. 

The key point of all wavelet methods is that the 
wavelet transform (WT) is able to discriminate structures 
as a function of the spatial scale, and thus is well suited 
to detect small scale structures on an image embedded 
within larger scale features. That is why, WT has been 
widely used for the structure analysis of galaxy clusters 
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(Slezak et al. 1994 Grebenev et al.||1995a[ [Rosati et"aL 



1995[ iBiviano et al.|[T996[ |Vikhlinin et al.||1997[ ). 

In the context of an individual IBIS/ISGRI observa- 
tion, we are interested in removing all large scale struc- 
tures from non-uniform sky background rather than in the 
detection of point sources. The task is greatly simplified 
by the fact, that the coded-mask aperture technique is 
not able to reconstruct an image of objects with a spatial 
size greater than the angular resolution of the telescope. 
For point source detection all structures more extended 
than a point source can be safely removed. In other words, 
we do not need any thresholds to discriminate noise and 
signal, we can remove systematic residuals with a given 
angular scale "as is" . In this way, the WT works as a non- 
parametric sky background approximation. In the similar 
way, the wavelet transform was used for subtraction of 
nonuniform background and for filtering images obtained 
with the coded-mask X-ray telescope ART-P aboard the 
GRAN AT mission ( jGrebenev et al.||1995bD . 



4.2. Wavelet decomposition method 

To decompose a sky image, we use the a trous digital 
wavelet transform (DWT) algorithm because it allows 



al. 


1994 


nel 


Kj 



easy reconstruction (Starck & Murtagh 1994 Slezak et 



Fj — Fj+i, where integer J is the so called 
spatial scale index. Each Fj is constructed by five weight- 
ing coefficients [1, 4, 6, 4, 1]/16 spaced by a 2'^~^ interval. 
Note, that each Fj can be roughly approximated by a 
Gaussian of width 2^~^ {Fi is a (5- function) . The convo- 
lution of an image with Fj preserves flux (HFj = l,Vj), 
and convolution with Kj emphasizes the structures with 
the characteristic size ~ 2^~^ pixels, or 2^+^ arcminutes 
in the case of the IBIS/ISGRI image (Ipix ^ 4'). Thus, 
low J values correspond to small spatial variations or high 
frequency, and high J reflects large spatial variations or 
low frequency. 

On the largest scale, N the kernel Kn = F^. The orig- 
inal image / can be easily decomposed to its convolutions, 
Wj ("wavelet planes of scale J"), with kernels Kj: 



I 



^j=iWj. (4) 

Therefore, we can consider Wj an image containing "flux" 
on scale J; the sum of all fluxes yields the original image. 
This is the basis of the wavelet decomposition algorithm. 
The original image is convolved with the wavelet of the 
first scale. The wavelet plane on the first scale is removed 
from the image. We then go to the next scale J + 1. In 
other words, we remove all small-scale features from the 
image before working at larger scales. That is why the 
small scale features (high spatial frequency) do not affect 
the convolution at larger scales (at low spatial frequency). 
Obviously, this algorithm greatly reduces the interference 
of the point source with large-scale kernels. 

In order to remove the large-scale systematic struc- 
tures seen in the left image of Fig. [S) we start with de- 
convolved source-free sky image of individual observation: 



Table 1. Reconstruction of the point source flux with 
different sets of wavelet planes. 

DWT scales, J Flux fraction Significance 



1 

1,2 

1-3 

1 -4 

1-5 

1-6 

1-7 



0.00435 
0.36121 
0.78648 
0.94208 
0.98549 
0.99647 
0.99913 



1.54 
127.95 
278.60 
333.72 
349.09 
352.98 
353.92 



original sky image 
1-8 



1.00000 



354.23 



wavelet scales used in this work 
1 -4,8 0.94210 



333.72 



1) flux image is decomposed to the wavelet planes. The 
systematic residuals are clearly seen at scales J = 5,6,7. 
After that, we 2) restore the original image by Eq. [4] omit- 
ting these scales, and 3) return the source fluxes to the 
sky as described in the general reconstruction method 
dKrivonos et~ar||2005[ [Goldwurm et al.||2003 ). 

We should stress, that DWT filtering is used on sky 
fiux images, where systematic residuals are emphasized. 
The sky variance map was not filtered, because it contains 
only formal uncertainties related to the exposure time for 
a given sky direction. 



4.3. Impact on point sources 

According to the above section, the cleaning procedure 
operates on the source-free images, i.e. images free from 
catalogued sources. Therefore, DWT filtering does not di- 
rectly affect the flux of known sources. 

However, if position of the source is unknown, its flux 
will not be removed from the detector before sky decon- 
volution, which means that source will appear on the sky. 
In this case, the DWT filtering will clean sky background 
with embedded point source. 

The important question is how the implemented DWT 
filtering affects the flux of such an point source. The main 
idea is that, the DWT procedure must not significantly 
change the point source flux. 

In order to investigate this issue, we used standard on- 
axis observation of the bright source Crab Nebula with 
total exposure 2.7 ks. We performed DWT decomposition 
of the reconstructed sky image of Crab, for the range of 
scales 1—8. The image was then reconstructed by summing 
up selected wavelet scales. Crab flux was estimated on the 
final sky image, convolved with the effective point spread 
function ( Krivonos et al.||2007bD . 

The Table [T] contains the measured Crab flux for the 
different sets of wavelet scales. The flux is expressed as 
a fraction relative to the flux in the original sky image. 
As seen from the table, ~ 94% of the point source flux 
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Fig. 8. Sky region around NGC 4151 accumulated with a 
sequence of 95 staring observations (spacecraft orbits 74- 
76). The total dead-time corrected exposure is 280 ks. The 
angular size of each image is 30° x 30°. The left mosaic 
image is obtained by the general method (Sect. |2|. The 
image on the right was produced by summing up images 
of individual observations corrected with a trous wavelet 
decomposition algorithm (Sect. [XT]). The standard devi- 
ation of source-free pixels on the left and right images 
relates as 1.3 and 1.0, respectively. The new hard X-ray 
source ICR Jl 1203+4531, detected during this observa- 
tion is labeled in green. To illustrate how the algorithm 
works, we extract the vertical profile from the green rect- 
angular region in the left image (Fig. [10|. 



resides in the high-frequency wavelet scales 1 — 4, and the 
rest contributes to the low- frequency scales 5 — 8. 

If we assume that the point source is significantly af- 
fected when 5% of its fiux is greater than the la survey 
detection threshold, then wavelet filtering distorts only 
the point sources at a detection level of > 20<j. Obviously, 
such strong known sources are removed from the detec- 
tor shadowgram before sky deconvolution, and returned 
to the sky image after the DWT filtering in steps 1-2. 
New sources detected in the survey, generally have a de- 
tection significance of less than lOa. We conclude, that the 
implemented wavelet filtering method does not affect the 
fiux of point sources, and does not introduce significant 
distortion to the survey sensitivity. 

4.4. Extragalactic sky 



We performed DWT cleaning procedure (Sect. 4.2) for 
each observation of the already mentioned NGC 4151 star- 
ing mode campaign. The resulting image is presented in 
Fig. [Sj right. It is clearly seen that all large-scale arti- 
facts are totally removed, leaving a clean sky image. In 
order to demonstrate the achieved improvement, we built 
a distribution of signal-to- noise ratios (SNR) for pixels on 
the source-free sky image. The SNR distribution of the 
original image (Fig. [sj left) is represented in Fig. |9] by 
the red histogram. It can be approximated by a Gaussian 
with cr ~ 1.3, which is consistent with the measured stan- 
dard deviation of image pixels. The SNR histogram of 
the wavelet filtered sky image is plotted in blue. The last 
is well approximated by normal distribution with a unit 



100.00 r 



10.00 



0.10 



0.0 




Fig. 9. Signal-to- noise ratio distribution of a number of 
pixels in a 280 ks staring mode deep observation of 
NGC 4151. The SNR distribution of the image obtained 
with the general method (Fig. [sj left) is shown in red. The 
blue histogram shows pixel variance distribution of the sky 



mosaic accumulated with WT cleaned images (Sect. 4.1). 



The dashed line represents the normal distribution with 
unit variance and zero mean. Assuming normal distribu- 
tion of = 183606 pixels on a cleaned image, we expect 5 
occurrences at a significance level a = ±4.2. However, we 
detect 12 and 2 excesses at negative and positive values, 
respectively. 



variance and zero mean, which means that the systematic 
noise has been completely removed. Thus, we conclude 
that, for the simplified case of extragalactic observation, 
the DWT filtering significantly supresses (practically re- 
moves) systematic noise. 

Till now, about 3 Ms of INTEGRAL exposure time 
has been performed in staring mode. Usually, these ob- 
servations are excluded from sky mosaics due to high sys- 
tematics. Cleaning these observations with DWT we can 
add them to the survey. 

It is interesting to note, that by averaging the 
archival staring observation of NGC 4151 we detected 
a new transient source IGR Jl 1203+4531 at sky po- 
sition R.A.=llh20m21.60s, Decl.=45d31m48.0s (equinox 
2000.0, uncertainty 4 arcmin). The source was found at 
the FOV edge where strong systematics prevented its de- 
tection before (see Fig. [s] and 10). The source is seen at 
S/N = 5.3 on the original mosaic with RMS = 1.3, 
which gives 4(j excess. On the cleaned sky the system- 
atics is gone, and the source has S/N = 5.7 on the 
image with RMS = 1.0, i.e. the source is revealed at 
5.7(7. This is demonstrated on the average image profiles 



shown in Fig. 10 The follow-up Swift/XRT observation 
of IGR Jl 1203+4531 revealed two nearby sources with 
coordinates R.A.=llh20m26.92s, Decl.=+45:34:53.77 and 
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R.A =llh20m33.76s, Decl.= +45:28:17.92 (error radii ac- 
cording to the "xrtcentroid" program are 5.96 and 5.31 
arcsec, respectively). 

4.5. Galactic Center region 

However, we are mainly interested in improving sensitiv- 
ity in the region of the Galactic Plane where most of the 
exposures were collected. The sky image of the GP with 
the maximum available exposure 20 Ms in GC) pro- 
duced by the general reconstruction method is shown on 
the upper image in Fig. 11 The sky background behind 



the bright sources is contaminated by strong systematics 
which significantly limits sensitivity for source detection. 
In the same data set, by taking into account the Galactic 
X-ray Background and using DWT sky filtering we ob- 
tain a new deep image of the GP which is demonstrated 
in the lower panel. As seen from the sky image, most of 
the systematic artifacts are removed, leaving a more or 
less uniform sky background. Obviously, the quality of the 
reconstructed sky is improved. 

To demonstrate the efficiency of the improved recon- 
struction method we built the SNR distribution of the 
source-free 30° x 30° region around the Galactic Center, 



shown in Fig. 13 The red and blue histograms represent 
general and DWT corrected sky image backgrounds, re- 
spectively. The SNR distribution in the former case has 
wide no n- Gaussian sidelobes. The blue histogram repre- 
senting the cleaned sky is narrow, but still far from normal 
distribution. This means that systematic artifacts are re- 
duced, but still present on the sky. We measured the stan- 
dard deviation of image pixels masking out bright sources. 
For the cleaned and general sky standard deviation re- 
lates as 1.33/1.84, which gives us the total sensitivity im- 
provement of ~ 28%. Taking out the irreducible Poisson 
statistics having unit standard deviation, the suppression 
of systematic noise in the ~ 20 Ms deep field of Galactic 
Center region is ~ 44%. 

5. Survey 

For our analysis we used all data publically available in 
July 2009 and observations performed as a part of the 
GRXE study program (PI Sunyaev). The latter is mainly 
based on the Russian quota of INTEGRAL observing 
time. 

The sky image of any individual observation was pro- 
duced by the modified sky reconstruction method de- 
scribed in this work. The obtained sky images were added 
to the all-sky mosaics covering the whole sky. 

The survey coverage area is shown in Fig. [12] We cal- 
culated a fraction of the sky area covered by the survey 
at nominal and effective sensitivity. In the first case, the 
sensitivity is essentially the detection threshold, estimated 
from actual exposure (black curve). The effective sensitiv- 
ity was estimated by multiplying the nominal error map 
by variances of local background. The general (Sect. |2]) 
and modified sky reconstruction method (Sect. [4|, are 
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^ 0.1000 



o 0.0100 



^ 0.0010 - 
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le-11 le-10 
5(j Limiting flux, 17 — 60 keV (erg s^^ cm^^) 



Fig. 12. Fraction of the sky surveyed as a function of the 
limiting fiux for source detection with 5<j significance. The 
black curve demonstrates sky coverage for nominal sensi- 
tivity. The effective sensitivity estimated for general and 
modified sky reconstruction methods are shown by the red 
and blue curve, respectively (see Sect.lsl). 
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Fig. 13. Distribution of formal pixel significance in a 30° x 
30° region around the Galactic Center (Fig. 11). For plot 
description see Fig. |9] 



represented by the red and blue curves, respectively. The 
limiting fiux and sensitivity for 10% and 90% of the sky 
coverage are summarized in Table [2] 

General sky reconstruction method: The survey's lim- 
iting fiux for the longest exposure is ~ 35% higher than 
nominal. For a large coverage area at high limiting fiux, 
which is typical for high-latitude observations, the effec- 
tive sensitivity is ~ 15% lower than expected. 

Modified sky reconstruction method: The survey effec- 
tive coverage is apparently closer to the nominal. In the re- 
gion of low limiting flux, the effective sensitivity is ^ 27% 
worse than nominal. The sky coverage at high flux is con- 
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Fig. 10. Illustration of a non-parametric approximation of sky background (a and b) and source detection (c and d) 
with an d trous DWT algorithm. The black histogram in the spanel (a) is a vertically averaged profile of the sky 
image in Fig. [Sj extracted from the denoted green region. Blue, green and magenta curves are vertical profiles of 
corresponding wavelet components at the spatial scales J = 5,6, 7 (~ 1,2,4 degrees). The sum of these components 
(red curve) represents sky image background approximation which is to be subtructed. The corresponding residual 
profile is shown in the bottom panel (b). The abscissa axis is measured in the offset from the position of newly detected 
hard X-ray source IGR Jl 1203+4531 (red arrow). The right panels (c and d) are made in the same way. The black 
histogram in the top panel (c) demonstrates an actual S/N pixel values extracted along the line perpendicular to the 
green rectangular region in Fig. [S) and crossing the source IGR Jl 1203+4531. The long dashed lines represent the 
actual detection thresholds in standard deviations scaled from the RMS value measured on all the image pixels 



sistent with that expected from an actual exposure, i.e. 
Poisson statistics. 

In Fig{T4]we show survey sensitivity as a function of ex- 
posure time. Generally, the survey sensitivity grows with 
exposure by T~^l^ as expected. The region of high expo- 
sure in the Galactic Center is noticeable by its reduced 
sensitivity (red dots) in contrast with extragalactic obser- 
vations (green dots) where the limiting flux is expected to 
follow pure statistics (blue dots). Again, the survey sensi- 
tivity of the modified method is closer to a nominal sen- 
sitivity. 



Table 2. Survey limiting flux in mCrab^ for source detec- 
tion at the 5<j significance level (17-60 keV, 1 mCrab = 
1.43 X 10"^^erg s"^ cm"^). 



Category Nominal Gen. method Mod. method 

Lim. flux 026 035 033 

10% sky 0.60 0.70 0.63 

90% sky 4.32 5.00 4.32 



6. Conclusion 

We presented the improved sky reconstruction method 
for the IBIS telescope which suppresses the systematic 
effects. Firstly, the method takes into account extended 
Galactic X-ray Ridge emission which strongly affects the 
background illumination of the ISGRI detector. Secondly, 
we applied a non-parametric (model-free) background ap- 



proximation based on an d trous wavelet decomposition. 
The wavelet cleaning was naturally integrated into the sky 
reconstruction process with the main advantage that we 
knew exactly what we were filtering out without distorting 
the original sky flux. The overall systematic noise in the 
~ 20 Ms deep field of the Galactic Center was reduced by 
~ 44% improving the the total sensitivity of observations 
by ~ 28%. The reconstructed sky images of high galactic 
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Fig. 11. Map of the sky region near the Galactic plane obtained with IBIS/ISGRI in the 17 — 60 keV energy band. The 
total exposure is about 20 Ms in the region of the Galactic Center. Upper panel: sky mosaic acquired by the general 
sky reconstruction method (see Sect.[2|. Lower panel: sky mosaic produced by an improved reconstruction algorithm 
(Sec t. [4| ). The corresponding signal-to- noise distributions of pixels in a 30° x 30° region around the GC are shown in 
Fig.fW 



latitude fields were practically free from systematic resid- 
uals and the sensitivity was consistent with that expected 
from Poisson statistics. 

Most of the INTEGRAL observing time was spent 
in the Galactic Plane and Center, giving us the possi- 
bility to conduct the most sensitive survey ever made of 
the Milky Way above 20 keV. The minimal detectable 
fiux with a 5cr detection level reached the level of 3.7 x 
10~^^ erg s~^ cm~^, which is ~ 0.26 mCrab in the 17- 
60 keV energy band. The survey covered 90% of the sky 
down to the flux limit of 6.2 x 10~^^ erg s~^ cm~^ (~ 
4.32 mCrab) and 10% of the sky area down to the flux 
limit of 8.6 x 10"^^ erg s"^ cm"^ (~ 0.60 mCrab). A cat- 
alogue of sources detected in the survey is presented in 
paper by Krivonos et al., (2010b, in prep.). 
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Fig. 14. Measured 5<j limiting flux for catalogued source 
positions on the final total average all-sky map, as a 
function of effective dead- time corrected exposure time. 
The top and bottom panels demonstrate the sensitivity of 
general and modified sky reconstruction methods, respec- 
tively. The blue points reflect nominal sensitivity chang- 
ing with exposure by T~^/^ as expected. The solid line 
represents the fit to the nominal sensitivity versus expo- 
sure: /f4 = 0.77 X (T/Ms)-^/2 mCrab. The green and 
red points are 5a actual sensitivity taking into account 
background variations in the field of each 20° x 20° sky 
projection (Sect.[5|. The sensitivity measured in (and out 
of) sky region |/| < 20° and \b\ < 15° is shown in red 
(green). 
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